Cavity cooling of a trapped atom using Electromagnetically-Induced Transparency 
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A cooling scheme for trapped atoms is proposed, which combines cavity-enhanced scattering and 
electromagnetically induced transparency. The cooling dynamics exploits a three-photon resonance, 
which combines laser and cavity excitations. It is shown that relatively fast ground-state cooling 
can be achieved in the Lamb-Dicke regime and for large cooperativity. Efficient ground-state cooling 
is found for parameters of ongoing experiments. 
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I. INTRODUCTION 

Control of the quantum dynamics of physical ob- 
jects is a prerequisite for quantum technological applica- 
tions. One important requirement is high-fidelity quan- 
tum state preparation. This often relies on cooling the 
physical system of interest to sufficiently low tempera- 
tures, such that a single quantum state, the ground state, 
is occupied with probability approaching unity. Laser 
cooling constitutes in this respect a successful technique, 
which is routinely used in the preparation stage of ex- 
periments with atoms and ions 0, Furthermore, the 
various laser-cooling concepts and schemes which have 
been proposed and partly tested over the years are be- 
ing considered for cooling molecules Q and more com- 
plex objects, such as for instance phononic modes of 
micromechanical resonators to ultralow temperatures Q 
and condensed-phase systems 0, @ . 

The idea at the basis of laser cooling of trapped par- 
ticles is to enhance scattering processes leading to a 
net transfer of the oscillator mechanical energy to the 
electromagnetic-field modes. Ground-state cooling of a 
harmonically-trapped particle, in particular, relics on a 
strong enhancement of the scattering processes which 
cool the oscillator, over the ones which heat the mo- 
tion. Several ground-state cooling techniques have been 
discussed that apply this basic concept in various ways 
[2, 0~E1 • One scheme that is relevant to the study per- 
formed in this work is known in the literature as EIT- 
cooling, where EIT stands for Electromagnetically In- 
duced Transparency jToj , and it uses coherent population 
trapping [lOL [Tl| due to quantum interference between 
laser-driven electronic transitions, in order to tailor the 
scattering cross section of the atoms to pursue ground- 
state cooling [T2], . The EIT-cooling scheme extends 
the basic concepts of velocity-selective coherent popula- 
tion trapping for free atoms [13, EH to trapped atoms, 
and has been experimentally demonstrated in Refs . flfjl - 
[l8j . Further studies can be also found in Refs. [1107 



' marc.bienert@physik.uni-saarland.de 



A further tool that is being increasingly considered in 
order to enhance scattering processes leading to trap- 
ping and cooling of atoms, is the strong coupling of one 
or more atomic transitions to a high-fincssc optical res- 
onator [22| - [30| . For such a system it has been predicted 
that quantum interference effects can e merg e from the 
quantized nature of the cavity mode [U [Hj. When in 
addition the field is interfaced with the atomic vibration 
in a trap via the mechanical action of light, further inter- 
ference effects can emerge that can increase the cooling 
efficiency H, 33 . The most recent observation of cavity- 
induced EIT I35l . l36j . and of mechanical effects associated 
with it [37l |38| . leads to the natural question whether 
and how EIT and cavity quantum electrodynamics can 
concur together to provide novel tools of control on the 
mechanical effects of atom-photon interactions. 

In this paper we present a theoretical study of the me- 
chanical effects of light on a trapped atom in a setup 
which supports cavity-induced EIT. We show that the 
combination of EIT and cavity quantum electrodynam- 
ics can give rise to a cooling mechanism, which for an 
accurate choice of parameters allows one to prepare the 
atoms in the ground state of the potential with probabil- 
ity approaching unity. Remarkably, high efficiencies are 
found for the parameters of the experimental setup re- 
ported in Ref. H3] ■ We show that the cooling dynamics 
can be often explained by means of a three-photon res- 
onance [HHH, which involves cavity and laser photons. 
For certain parameter choices cooling results from inter- 
ference in the mechanical effects of the atomic interaction 
with the cavity and laser fields. 

This article is organized as follows. In Sec [IT] the physi- 
cal system and the theoretical model are introduced. The 
assumptions are discussed, which are at the basis of the 
theoretical treatment in this article. In Sec. IHII the dark 
resonances and dressed states of the Hamiltonian for the 
electronic and cavity levels are reported, and the cavity 
and atom excitation spectra are discussed. In Sec. IIVI 
the basic equations of cooling are derived, and in Sec. 
IVl the corresponding predictions are reported for exper- 
imental parameters based on Ref. [37} • The conclusions 
are drawn in Sec. IVI| while in the appendices the details 
of the calculations in Sees. HTl IlIII and IIVI are presented. 
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FIG. 1. a) Setup of the system. An atom is confined inside an 
optical resonator by a harmonic trap with frequency v. The 
atom is transversally driven by a laser at Rabi frequency Ql 
and couples to the cavity mode with vacuum Rabi frequency 
g. The cavity mode is pumped by a laser (coupling strength 
Qp) and decays with rate k. (4>c) denotes the angle be- 
tween the axis of the motion and the laser (cavity) wave vec- 
tors, b) Relevant electronic transitions. The transverse laser 
(cavity mode) drives the transition \gi) — » |e) (|g2) — ► |e)). 
The excited state \e) decays spontaneously into the stable 
states and \g2) with rates 71 and 72, respectively. 



II. THE THEORETICAL MODEL 

A single atom of mass M is confined inside of an op- 
tical resonator by an external harmonic trapping poten- 
tial and is irradiated by a laser field, while the cavity 
is pumped by a second laser field at strength fip. Two 
atomic dipolar transitions couple to the laser and the 
cavity mode, respectively, and share the same excited 
state, forming a A-shaped configuration of levels. The 
atomic center-of-mass motion is treated in one dimension 
along the x axis. Although the one-dimensional treat- 
ment seems to be a strong restriction, in the Lamb-Dickc 
regime which we will assume here, the rate equations for 
the cooling dynamics can be split into three independent 
sets of rate equation, one for each direction of motion. 
Thus, in this limit each spatial dimension can be treated 
separately [43| . Figure Q] illustrates the setup and the 
level configuration, highlighting the geometry of the laser 
beam and the cavity axis with respect to the axis of the 
motion. 

In the following we report some of the relevant param- 
eters and introduce the basic notation. The atomic level 
configuration is reported in Fig. QTb). The atomic levels 
are denoted by the stable states |<?i) and ^2} which are 
coupled by a dipolar transition with moments pi and p2 , 
respectively, to the common excited state |e). We denote 
the level frequencies by ui\, u>2 & n d w e . The atomic dipolc 
transition |<7i) O |e) is driven by a laser with frequency 



wl and Rabi frequency f^L- The transition 1 52) <-> |e) is 
coupled to a mode of the optical resonator with frequency 
wc, lincwidth 2k and vacuum Rabi frequency g. The ex- 
cited state has radiative linewidth 7 and decays into the 
state \gj) with rate jj (j = 1, 2) such that 7 = 71+72 [H|. 

The atomic center-of-mass motion is confined by a har- 
monic potential of frequency v which is independent on 
the electronic state of the atoms. This situation can be 
realized when the single particle is an ion in a Paul or 
Penning trap p5| . or a neutral atom when the confin- 
ing potential is a dipole trap under specific conditions 
(such as the magic wavelength) [4(| |47| . The geometry 
of the setup is fixed by the angles </>l and 0c which give 
the orientations of the laser and cavity wave vectors with 
respect to the axis of motion, see Fig. [TJ 

The dynamics of the system include the mechanical 
coupling of the atom with the electronic transition via 
absorption and/or emission of photons. Moreover, it de- 
scribes the coupling of the cavity mode with the lon- 
gitudinal modes of the electromagnetic field by the fi- 
nite transmittance of the mirrors. Correspondingly, the 
Hamiltonian H. governs the dynamics in the Hilbert space 
of the degrees of freedom of the atom, the cavity mode, 
and external modes of the electromagnetic field. It can 
be decomposed into the sum 

% = H + W 1 + W K + ffemf • (1) 

Here, H is the Hamiltonian for the dynamics of the sys- 
tem, composed by cavity mode, electronic bound states 
and center-of-mass motion of the atom, and contains 
the coupling with the lasers, which are treated as classi- 
cal fields. Term H om { describes the Hamiltonian of the 
modes of the electromagnetic field external to the res- 
onator. The coupling between the atomic dipole transi- 
tions and these modes is given by W 7 , which is respon- 
sible for the radiative instability of excited state |e) at 
rate 7. The longitudinal modes of the external electro- 
magnetic field couple to the cavity mode via the finite 
mirror transmittance. This coupling is incorporated by 
the term W K and gives rise to cavity losses at rate k. In 
what follows we discuss in detail the form of the individ- 
ual terms. 



A. Hamiltonian of the atom-cavity system 

The Hamiltonian 

H = H cxt + H int + H c&v + W (2) 

governs the dynamics of the composite system of atom 
and cavity. In a frame rotating with the lasers' fre- 
quencies, the individual terms on the right-hand-side of 
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Eq. © read 



(3) 



-Hint = -M c2 |e)(e| + h[8i - S c2 ]\gi){gi\ + hA\g 2 ){g 2 \, 

(4) 

H cav = — ftAa^a, (5) 

and describe the center-of-mass motion of the atom in 
the harmonic potential of frequency v, the unperturbed 
dynamics of the internal electronic states, and the dy- 
namics of the cavity mode, respectively. The operators b 
and W in Eq. © annihilate or create single vibrational 
excitations and are connected with the atom's position 
operator by the relation 



(6) 



where £ = \J h/2Mv denotes the size of the ground-state 
wave packet. The operators a and a) in Eq. (JSJ) annihi- 
late and create, respectively, a photon in the considered 
mode, and 



A = up — luq 



is the detuning between the cavity and the probe. The 
other detunings occurring in Eq. ^ are 



*i = 

The interaction part 



CJ L - (w e - UJi) , 
WC - (w e - L0 2 ). 



W = W P + Wl + W c 



is composed of 



Wp = — - \a + a'] 



W h (x) 



[\e){ gi \e lkxcos ^ +H.c] 



W c (x) = hg{x)[\e)(g 2 \a + R.c. 



(8) 
(9) 



(10) 



(11) 

(12) 
(13) 



where Wp is the drive of the cavi ty field b y the probe of 
power P with strength Qp = 2^ Pn/hup, Wl{x) is the 
coupling of the atomic dipole \g\) -o- |e) to the control 
laser with Rabi frequency and Wc(x) the Jaynes- 
Cummings interaction between the cavity mode and the 
transition \g 2 ) <-> e) with coupling constant 

g{x) = g cos(kx cos 0c + ^p) ■ 

The quantity k is the wavenumber of the control laser 
field and of the cavity mode [Hj]. Moreover, tp determines 
the equilibrium position of the atom in the trapping po- 
tential with respect to the cavity mode function. 



B. Coupling to the electromagnetic field external 
to the resonator 



We include the electromagnetic field outside the cavity 
into the theoretical model using a Hamiltonian descrip- 
tion instead of a master equation. This provides more 
insight into the basic processes underlying the cooling 
dynamics. The Hamiltonian 



EC72) 



E(«). r 



(14a) 



accounts for the energy of the transversal and longitudi- 
nal field modes which couple independently to the atomic 
dipoles and the cavity, respectively. The superscript ~fj 
(k) indicate that the sum is restricted to the modes which 
couple quasi- resonantly with the transition \gj) — > \e) 
(with the cavity mode) and the Hamiltonian is reported 
in the reference frame at which the corresponding tran- 
sition rotates at the frequency of the driving laser. The 
quantity lo^ = c\k\, with the speed of light in vacuum c, 
denotes the frequency of the electromagnetic field mode, 
and cj e is the annihilation operator of the mode la- 
beled by the wavevector k and polarization eg _L k with 
e = 1,2. 

The interaction between the atomic dipole \gj) f-> |e) 
and the external electromagnetic field reads 



W 7 (x) = W~ /1 (x) + W j2 (x) 



3 =1 . 2 fe.e 



g { y}\e){gM {% - S * )x ct e + K, 



(15) 



where the coupling constant ffjj = Pj ' is pro- 

portional to the scalar product of the atomic dipole mo- 
ment pj with the vacuum electric- field £^e^ . The cou- 
pling of the cavity mode to the longitudinal modes of the 
external electromagnetic field is given by 



W K 



\ ( K ) t 



H.c 



(16) 



with coupling constant g^ c (a detailed form in terms 
of the physical parameters can be found for instance in 
Ref . [4!| ) ■ We remark that the operator W K does not de- 
pend on the atomic motion, but describes dissipation pro- 
cesses which will be instrumental for cooling the atomic 
motion [24[. We furthermore emphasize that quantum 
noise due to fluctuations of the atomic dipoles and the 
cavity field is systematically incorporated in the theoret- 
ical description by the couplings Eqs. (fT5j) and (fTBT> . 
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C. Basic assumptions and perturbative expansion 

Throughout the paper we focus on the regime where: 
(i) The atom is tightly confined: The size of the atomic 
wavepacket Ax is much smaller than the lasers' wave- 
lengths. This is the so-called Lamb-Dicke regime pi. [50j. 
It is found when t he inequal ity kAx <C 1 is satisfied. Us- 
ing that Ax = £ t ^j2(m) + 1 by averaging over a thermal 
state with mean occupation number (m), the inequality 
can be rewritten as 



77v/2(m) + 1 < 1 



(17) 



where 77 = fc^ is the Lamb-Dicke parameter, which scales 
the mechanical effects of light on the atomic motion [5(| ■ 
For later convenience we define 

t?l = V cos L , 
77c = r\ cos 4> c , 

which account for the geometry of the mechanical cou- 
pling of laser and cavity to the axis of motion, respec- 
tively. 

We further assume that (ii) the laser driving the cavity 
is sufficiently weak such that the average photon number 
of the cavity mode is much smaller than unity. This 
corresponds to taking the small parameter 



«1, 



(18) 



which gives the mean number of intracavity photons 
when no atom is present. The requirement (fT8|) can be 
achieved for high-finesse cavities by adjusting the param- 
eters f2p and A accordingly. 

This regime allows for a perturbative treatment of the 
dynamics. We expand the total Hamiltonian into differ- 
ent orders of the Lamb-Dicke parameter by performing 
a Taylor expansion in rj of the exponentials cxp[. . . x] in 
Eqs. (H2J) and (TJS]) and the function g(x) in Eq. {JT3J) (Hj. 
Moreover, we separate the weak driving term Wp, i.e., 
the first order term in |e|, from Ho, so that the Hamilto- 
nian takes on the form 

H = H + W P +H 1 +O{if). (19) 

The term Ho is in lowest order in 77 and |e|, and reads 

Ho = H + W 7 (0) + W K . (20) 

Here, 

Ha = H cxt + H opt + -f/cmf , (21) 

and 



H, 



O [ ) t 



W L (0) + W C (0), (22) 



which is given in zero order in the perturbative expan- 
sion. In this order of perturbation theory, the stable state 
of the atom-cavity system is 



|*Bt> = |ff2,0> 



(23) 



which is the product state of the atom in \g%) and no cav- 
ity photon, \n = 0). Under the influence of H (assuming 
the regime in which the Wigner-Wcisskopf approximation 
can be applied [49j). the atom-cavity system relaxes into 
the state l^st), while the atomic center-of-mass motion 
evolves coherently and is decoupled from the electronic 
dynamics. 

Coupling between the center-of-mass motion and the 
light is introduced in first order in the Lamb-Dicke pa- 
rameter rj. The corresponding term of the total Hamil- 
tonian reads 



Hi=Fx+ F 7r 
3=1,2 



(24) 



It accounts for the mechanical effects due to absorption 
and emission of photons from/to the control laser and 
the cavity mode, where the operator F 



F =F L + F c 



-- (W h (x) + Wc(x))\ x=Q ., 



can be interpreted as a force operator [52j, while 



Fr 



VW.Ax) 



x=0 



(25) 



(26) 



gives rise to the stochastic force, which is associated with 
the recoil due to spontaneous emission J52j. Here, e x is 
the unit vector along the axis of the motion. 

The Lamb-Dicke regime is found when the condition 
set by Eq. (fl~7|) is fulfilled, and it allows one to assume 
the separation of time scales which determine the center- 
of-mass motion of the atom, and the common, internal- 
state dynamics of the atom-cavity system. At lowest or- 
der in the Lamb-Dicke expansion, the internal and exter- 
nal atomic degree of freedom evolve independently The 
relevant time scale is given by the smallest relaxation 
time of the atom-cavity system. At higher order, the 
atomic center-of-mass motion experiences the force due 
to the gradient of the electromagnetic field over the ex- 
tension of its wave packet. These processes take place on 
a time scale, which is slower by a factor ry 2 with respect 
to the time scale of the internal motion. In this regime 
the internal motion follows the external motion adiabati- 
cally [ill [H| . It is therefore instructive to first study 
the internal dynamics, neglecting the coupling with the 
external motion: This permits one to identify scattering 
processes which leads to cooling. 



III. LEVEL STRUCTURE AND PROPERTIES 
OF THE ATOM-CAVITY SYSTEM 

We now discuss the dressed states of Hamiltonian H opt , 
Eq. (|22|) . Condition (fT8|) allows one to restrict the cavity- 
atom Hilbert space to states which contain at most one 
excitation of the cavity mode. The dynamics then takes 
place within the subspacc spanned by the states 

{| 52 ,l),| e ,0),| gi ,0),| 52 ,0)}. 



5 



These coupled levels constitute an effective four-level sys- 
tem, which is depicted in Fig. [2l Such configuration of 
levels has been studied in the literature for the case in 
which the states are four different electronic levels, that 
are coupled by classical laser fields (39144^ . |54| . For this 
case it has been found that this level scheme can exhibit 
dark resonances due to quantum interference between ex- 
citation paths involving three photons. The purpose of 
this section is to analyze the spectroscopic properties of 
the four-level system arising from the atom-cavity cou- 
pling, in the regime in which we can neglect the coupling 
with the center-of-mass motion. We will focus, in par- 
ticular, on the conditions under which dark states exist 
as they will turn out to be relevant for the cooling dy- 
namics, and refer the reader to Ref. |42j for an extensive 
analysis. 

The eigcnspcctrum of H opt , Eq. (1221) . within the con- 
sidered subspace is displayed in Fig. [3] as a function of 
the frequency of the transverse laser, under the assump- 
tion that the state |<7 2j 0) is weakly coupled to the other 
states. The frequencies of the states e, 0), \gi, 0), \g2, 1), 
which diagonalize H opt when the coupling to laser and 
cavity is set to zero, are indicated by the dashed curves. 
In presence of the coupling with the fields, the frequen- 
cies are shifted and degeneracies become avoided cross- 
ings. The resulting eigenfrequencies are u>±, lu , which 
correspond to the eigenvectors |±) and |o), being su- 
perposition of the three states |e,0), |<7i,0), \g 2 ,l)- The 
horizontal part of the curves, in particular, denote the 
frequency of the dressed states of the Jaynes-Cummings 
Hamiltonian, namely superposition of states \g2, 1) and 
e, 0). These states appreciably mix with state |<7i,0) at 
the level crossing with the curve 5\ — 5 C 2- This is also vis- 
ible by inspecting the linewidth, indicated in the figure 
by the breadth of the curves. 

We now analyze the excitation spectra of cavity and 
atom, namely, the rate of photon emission by the cav- 
ity and the atom, as a function of the probe frequency 
A. They are evaluated for the stationary state of the 
(internal) atom-cavity system, and are proportional to 
the probability that the cavity contains one photon and 
that the atom is in the excited state, respectively. In Ap- 
pendix A a detailed calculation is reported, which shows 
that they take on the form 



a) 



5 c K xc (A)cx|J- K (A)| 2 , 
S£ c (A)oc |J- 7 (A)| 2 



(27a) 
(27b) 



with 



7- K (A) = ^ + A -^^ + A + ^]-f ! 



/(A) 



J- 7 (A) = g- 



A - Si 



/(A) 



(28b) 



where the superscript k (7) indicates that it refers to 
cavity (atom) emission, see App. [A] for details. The rel- 
evant resonances can be identified with the poles of the 
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FIG. 2. Level scheme of the atom-cavity system in a) lab- 
oratory and b) rotating frame. The vertical axis gives the 
frequency w. Downward arrows denote negative detuning. In 
b) the dressed states |<?i,±) and the states |<?2,e') with their 
corresponding frequency shifts are reported (see text). The 
detunings are defined in Eqs. ©-(O) and ((29)) . 



function /(A), that correspond to the frequencies of the 
dressed states of H opt . 

The excitation spectra for cavity and atom are dis- 
played in Fig. 2k) an d b) , respectively, as a function of 
A. We recall that in absence of the atom, the cavity 
excitation spectrum is a Lorentz curve centered at the 
mode frequency with full width 2k, and corresponds to 
the grey curve reported in the figure. The solid curve cor- 
responds to the cavity excitation spectrum in presence of 
the atom. The three peaks are identified with the dressed 
states of the system. Their centers are marked by the ar- 
rows in the abscissa and correspondingly by the arrows in 
Fig. [3] These peaks are also present in the atomic exci- 
tation spectrum, Fig. UJd). Additionally, we observe that 
the spectrum vanishes at a value of the probe frequency, 
and in this frequency region it exhibits a Fano-like profile. 
This is verified when the three-photon resonance condi- 
tion is fulfilled, namely, when states \g 2 , 0) and \g\, 0) are 
resonantly coupled, and corresponds to 



3 TP 



A - 61 = 0. 



(29) 



Moreover, the cavity excitation spectrum exhibits two 
points, where the cavity response is minimal. These min- 
ima are marked with two circles in Fig. Hk). 

The Fano-like profile and the minima in the cavity exci- 
tation spectrum can be understood in terms of quantum 
interference between the dressed states and state \g2,0), 
which can give rise to significant effects even though state 
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FIG. 3. Frequencies and linewidths of the dressed states of 
the atom-cavity system as a function of Si. The frequencies 
are found diagonalizing H opt , Eq. ((22} , over the basis of the 
unperturbed states {|e,0}, |<?i,0), \g2, 1}}, and correspond to 
the centers of the black curves. The corresponding linewidths 
determine the curves' width. The horizontal dashed-dotted 
lines and the dashed diagonal line give the frequencies of the 
unperturbed states as a function of Si . The arrows mark the 
frequencies of the resonances shown in Fig. [4] The parameters 
are k = 2u, 7 = IO17 g = 20u, = 12^, tp — 0, <5 C 2 = 20v. 

I <72,0) is very weakly coupled to the other levels. In or- 
der to gain more insight, we first analyze the atom cavity 
levels in a unitarily equivalent scheme, where H cav + Wp 
is diagonal. We consider the states 

\ 9j ,e')=D c (e')\ gj ,0) (j = l,2) (30) 

with the displacement operator D c (e) = cxp(ea^ — e*a — 
iA\e\ 2 t) for the cavity mode and e' = f2p/2A. For 6 c2 + 
A — 5i = 0, the states Eq. (f3TJ|) are resonantly coupled, 
such that 

\D)^e'gcos<p\g u e')~^\g 2l e') (31) 

is an eigenstate of the Hamiltonian H opt , that has zero 
projection onto the electronic excited state. In the limit 
in which the pump is weak and out of resonance the state 
\D) is to good approximation the stationary state of the 
system and is a dark state. Equation (|29|) then gives 
the correct resonance condition and explains the posi- 
tion of the minimum in the atomic excitation spectrum. 
Similarly, the cavity excitation spectrum exhibits min- 
ima when the state 1 1/ 2) 0} is resonantly coupled to one 
of the dressed states \gi,±), which diagonalize the laser 
interaction H mt + Wl(0). 

Another dark resonance is found when the states \g 2 , 1) 
and 0) are resonantly coupled, namely, when <5i = 5 C 2- 
We denote this situation by "two-photon resonance" . In 
this case the state 

\D M ) oc 5 cos ^|.gi,0) - -y 1^2,1) (32) 

is an eigenstate of the Hamiltonian H opt and is stable 
over a time scale, in which cavity decay can be neglected. 
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FIG. 4. Excitation spectra of a) cavity and b) atom, in 
arbitrary units, as a function of the probe frequency A, in 
units of v, for <5 C 2 = 20^ and Si — lQu (black solid line) and 
Si = S C 2 = 20v (dashed line). The gray curve in a) gives 
the Lorentzian excitation spectrum of the unperturbed cavity 
and is plotted for comparison. The spectra consist of three 
peaks at the frequencies of the dressed states of the manifold 
{|e, 0), \gi, 0), |<5)2, 1)}, which are marked with arrows. The 
circles mark the frequencies for which the excitation spec- 
tra exhibits local minima, corresponding to approximate dark 
resonances. In b) the atomic excitation spectrum vanishes at 
A = Si — S C 2 (small circle), due to three-photon resonance 
between |gri,0) and |(/2,0). The other parameters are k — 2v, 
7 = lQu, g = 2Qu, n L = 12^, (p = 0. 



This situation is reported by the dashed curves in Figs. 
|4]a) and b). In this case, the atomic excitation spectrum 
exhibits only two peaks, and the dark resonance is at the 
frequency where the cavity output is maximal. 

In the next section we take into account the motion 
of the atom, and trace back the properties of the rates 
of heating and cooling transitions to the characteristics 
identified in the excitation spectra. 



IV. THEORY OF COOLING IN THE 
LAMB-DICKE LIMIT 

In the Lamb-Dicke regime, the atom's external and in- 
ternal degrees of freedom are weakly coupled. In this 
limit, one can assume that the system, composed by 
cavity and electronic excitations of the atom, reaches a 
steady state over a time scale which is much faster than 
the one at which the motion evolves. The corresponding 
theoretical treatment has been discussed in extent, for 
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instance in Rcf. [13J, [53|, [55| . This justifies the formula- 
tion of the dynamics of the external degrees of freedom 
in the form of a rate equation for the occupations p m (t) 
of the vibrational state |m), which reads 



dp,, 



dt 



=(m + l)j4_p TO+ i - [(m + 1)^4+ + mi_]p„ 



m — 1 



(33) 



The average phonon number obeys the equation (m) (t) = 
—T(m)(t) + A + with the cooling rate 



T = A_ - A+ 



(34) 



When r > 0, i.e., > A + , a stationary state exists. At 
steady state, the flow of population fulfills the detailed 
balance condition, and the mean occupation at steady 
state reads H6| 



A 



with mean phonon occupation 



(m) s 



A 4 



A_ - A+ 



(35) 



(36) 



These formula show that high cooling rate are reached by 
maximizing the value of A_ as well as the ratio A—/A+, 
while the ratio A_/A + alone controls the final tempera- 
ture. Various cooling schemes resort to different control 
tools in order to achieve either fast cooling and/or low 
temperatures. In the following we will characterize the 
dynamics which can lead to efficient ground-state cooling 
in this setup. 



A. Cooling and heating rates 



The transition rates A± in Eq. ([53")) can be calculated 
in perturbation theory using the resolvent of the Hamil- 
tonian. We refer the reader to App. [A] for the detailed 
evaluation. They describe absorption of a photon of the 
probe laser pumping the cavity, followed by emission into 
the modes of the external field either by atomic or by 
cavity decay, and can be written as the sum of various 
contributions of different physical origin: 



±|2 



r=l,2 
ft 



(37) 



The quantity D is a diffusion coefficient, originating from 
the mechanical effects of the spontaneously emitted pho- 
ton, while Tp ,=t denote the transition amplitude associ- 
ated with the absorption (— ) or emission (+) of a vi- 
brational phonon by either scattering the probe photon 
into the external modes by cavity decay (£ = k) or by 
spontaneous emission along one of the two transitions 



(I = 71.2)- The subscript F indicates whether the me- 
chanical effect is due to the laser photon (F = L) or to 
the cavity photon (F = C). For a given scattering pro- 
cess, these mechanical effects can interfere constructively, 
as one observes from the fact that the corresponding con- 
tribution, |7^' + Tq^I 2 , is the coherent sum of the in- 
dividual amplitudes. A schematic representation of the 
corresponding scattering processes is depicted in Fig. [5] 
A similar decomposition, as the one in Eq. (|37j) . was 
identified for the heating and cooling rates of a trapped 
atom, whose two-level, dipole transition couples simul- 
taneously to the cavity mode and to a transversal laser 
[55[ | . In such a setup several interference processes have 
been identified which can be tuned in order to enhance 
the cooling efficiency. In the present system, where an 
additional electronic transition is involved, a further in- 
terference process can take place, which is analogous to 
EIT. 

Below we discuss in detail the individual terms for the 
system considered in this work. In the following we ne- 
glect the decay along transition \e) — > |<ji), and take 
72 = 7, unless otherwise specified. This assumption leads 
to a considerable simplification of the analytical expres- 
sions we are going to report and does not affect quali- 
tatively the cooling dynamics for the parameter regime 
here explored, as numerical checks show. 

We first consider the diffusion term D in Eq. (|37p . 
This term describes scattering processes where a probe 
photon is spontaneously emitted by the atom, and where 
the mechanical effect is due to photon recoil by sponta- 
neous decay. It can be written in the form 



^ = 7W 2 |Td| 



(38) 



where W2 is a geometrical factor which depends on the 
atomic dipole pattern of radiation and 



ftp 

7d = -i?7— cos(/3j" 7 (A) 



(39) 



is the transition matrix element (apart of some constant 
factors), where the function .F 7 (A) is given in Eq. (|28bjl . 
Hence, D is proportional to the atomic excitation spec- 
trum in Eq. (|27a|l . and thus to the probability that the 
excited state is occupied. In particular, it vanishes at 
a node of the cavity mode function where cos<p = 0, 
since there, the state I52) docs not couple to the ex- 
cited state in zero order in the Lamb-Dicke expansion. 
Moreover, for general values of ip, D becomes zero when 
the atomic excitation spectrum vanishes, which is verified 
when A = Ao, with 

A = 5i - 5 c2 , 

namely, when states \g2, 0} and \g±, 0} are resonantly cou- 
pled. For these parameters the population of the atomic 
excited state is zero because of destructive interference 
between the excitation paths. Consequently, diffusion 
can be suppressed. This property is also at the basis 
of the so-called EIT-cooling mechanism 13j. Differing 
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from EIT cooling, suppression of diffusion is here due to 
a three-photon interference process. 

We now turn our attention to the transition amplitudes 



V: 



fip fi^ v(in + A =p v)g cos ip 



2 4 

rip si? 



/(A T ^)/(A) 



vg 2 cos 2 <p 



2 4 /(A T z/)/(A) 



(40a) 
(40b) 



They describe the scattering processes where absorption 
or emission of a photon of the transverse laser leads to 
a change in the vibrational motion and are depicted in 
Fig. [SJd) and c). These scattering amplitudes vanish at 
the node of the mode function where cos <p = 0, and when 
the control laser is perpendicular to the axis of motion 
= 0), namely, when there is no mechanical effect 
of the control laser along the direction of the motion. 
When this is not verified, these amplitudes may become 
maximal when the parameters are chosen so that the en- 
ergy defect of the intermediate scattering states becomes 
minimal. Cooling transitions are enhanced over heating 
transitions, in particular, for the parameters for which 
equation 



B. Three-photon resonance: EIT cooling and 
beyond 

In this section we study the cooling dynamics under 
the assumption of three-photon resonance, 

<5 T p = A + S c2 - £i = , 

corresponding to setting A = Ao- In this case the transi- 
tion rates simplify considerably: Due to destructive inter- 
ference of the excitation paths, the excited state |e) is not 
occupied and the diffusion term D in Eq. (|37[) vanishes. 
This property also leads to a considerable simplification 
of the transition amplitudes T c ,=l= , since J-" 7 (Ao) = 0. In 
this case the rates take on the form 153 



A± =|e|Vff 2 7>< 



1 + C± 



£ (1 + C ± f + (°jL-u±S 1 + ^C ± (v T A) 



where 



T 



rjl cos 2 tp ■ 



Vc sin2 <P 



(43) 



(44) 



contains the dependency on the geometry of the setup, 
and the parameters 



Re{/(A + ^)} = 

is fulfilled. 

Finally, the transition amplitudes 



Tc'* = -f/c^-sin^J^AT^) x 



(41) 



[gcos 2 ipT 1 (A) + {iK + AT^)J r 4A)] , (42a) 
~ ± Q P sin2(£ 
T c = - 71c — ^ X 

, 9 [J- 7 (A =f v) T K {A) + J- 7 (A) J" K (A =f u)] (42b) 



describe the scattering processes where absorption or 
emission of a photon of the cavity field leads to a change 
in the vibrational motion. They are depicted in Fig. [5Jd) 
and c). These transition amplitudes vanish when the cav- 
ity axis is perpendicular to the axis of motion (77c = 0), 
namely, when there is no mechanical effect of the cavity 
photon along the direction of the motion, or alternatively 
when the atom is situated at an antinode of the mode 
function, so that the derivative of the cavity-ion poten- 
tial vanishes and the corresponding force operator is zero. 
Term T£ !± , Eq. (|42b| , also disappears at a node of the 
cavity standing wave. Interestingly, both terms are prod- 
ucts of the Fano-like factors T 1 , Eq. (|28b|) , connected 
with the atomic excitation spectrum, and J 7 ^, Eq. (|28ap . 
which reflects the characteristics of the excitation spec- 
trum 5 e K xc , Eq. (j27bl) . 

We now identify some conditions, for which ground- 
state cooling can be implemented. 



C± = C- 



(A =F v) 2 

are proportional to the single-atom cooperativity 



C 



g 2 cos 2 ip 
«7/2 



(45) 



(46) 



The rate of cooling (heating) is maximum, when state 
|<?2 ; 0) is set on resonance with the red (blue) sideband of 
one of the dressed states of H op t. Cooling is enhanced, for 
instance, when the denominator of rate A_ in Eq. (|43|) 
becomes minimal, which is satisfied when relation 



is fulfilled. 



v+(v + A)- 



g 2 cos 2 (p 



{A + is) 2 



(47) 



1. Cavity EIT cooling 

In order to acquire a better understanding on the role 
of the resonator and identify the regimes in which cooling 
occurs, we first analyze the case C± <C 1. The parameters 
C± become small for (i) small cooperativities C, (ii) when 
the pump on the cavity is far-off resonant, or (Hi) for 
k <€. |A± v\. Then, the rates are given by 



A<<=|e|Vx 



2 9 

19 v 



v 2 \ 



— - u{vtoi) + — ^ K i +{ATV 



(48) 
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FIG. 5. Schematic representation of the scattering processes leading to a change of the motional state, corresponding to 
Eqs. (|39[) . l|lQ| an d (|42|l . The arrows show the sequence of processes leading to a transition |^ at ,m) — >• |\t s t,m± 1). The 
intermediate internal states are denoted by the rounded box. Process in (a) terminates by a spontaneous decay, in which the 
mechanical effect is due to the spontaneously scattered photon. In (b) and (c) the transition is terminated by emission of 
a photon by spontaneous emission and cavity decay, respectively, and the mechanical effects are due either by absorption or 
emission of a photon of the control laser (F = L) or of the cavity mode (F = C). 



For asymptotically small values of C, case (i), or large 
dctunings A| ^> v, re, case (ii), the rates assume the form 
of 



|2 — 2 

e| rj - 



1 1 

19 v 



(49) 



of {m) s tr.min = (^/2^) 2 , which can be very small provided 
that that the cavity lincwidth is much smaller than the 
trap frequency. 



2. Cavity EIT cooling in the strong coupling regime 



whose functional dependence on the parameters is the 
same as in EIT cooling (l3j |. For analyzing case (in) we 
use the definition of C, Eq. (|46|) . yielding 



.4 



|e|V- 



2 9 

19 v 



v 2 % + 



(50) 

For A = it corresponds to EIT cooling with a modified 
Rabi-frequency f2 2 /4 -> + g 2 cos 2 ip. 

All cases here discussed correspond to a mean phonon 
number at steady state 



(m) s 



[Q 2 /4-^(^ + (5i)] 2 +7 2 ^ 2 /4 
4^i(0 2 /4-r/ 2 ) 



(51) 



where Q denotes the corresponding (modified) Rabi fre- 
quency. It achieves the minimum value, (m) s t,min = 
(7/4<5i) 2 , when 2 /4 = v(v + Si). The correspond- 
ing cooling rate at this value reads W = A- — A + ~ 
|e| 2 ^4 5 2 / 7 . 

Resonant driving, A = 0, also leads to EIT-like cool- 
ing, even for arbitrary values of the cooperativity C, what 
is especially interesting for Cq = C±| A=0 3> 1. Then, the 
transition rates between the vibrational levels are given 
by 



A\ = 



1 2 — 2 
e\ 77 - 



ig 2 v 2 



,2 7_' 



(52) 



with the modified Rabi-frequency ^ 2 /4 — > 



7V 2 /2re and the modified atomic linewidth 7' = 7C0. 
For asymptotically large Co this case leads to the mini- 
mal mean vibrational occupation number at steady state 



We now consider the case of large cooperativities, 
C» 1, and small cavity linewidths, such that re <C v. 
Moreover, we assume that the system is driven at three- 
photon resonance, 5tp = 0. 

When the pump is tuned on the blue sideband transi- 
tion of the cavity, A = is, one encounters the situation 
C+ = C and C- w C^, fulfilling C + > C_. The result- 
ing transition rates between the vibrational levels take 
on the form 



|e|V 



C 7 ' 



A_ 



1 2 — 2 



5 2 7(1 + C_) 



f (l+C_) 2 + (<5 1 -<5 opt ) 5 



(53) 
(54) 



where 



(5 opt = ^L/(4^)+ff 2 cos 2 <p/{2v) 



This corresponds to the case where the weak pump is in 
resonance with the red sideband of the dressed state |+) 
of -ffopt- It becomes clear from Eq. (|53|) that the cooling 
is larger than the heating rate by the factor 



A. 



C~ 



i,y- 



Additionally, diffusive processes are suppressed due to 
three-photon resonance. In this regime, the mean occu- 
pation number reads 



(m) st = 



L (l + C^) 2 +4(,5 1 - ( 5 opt ) 2 /7 2 



, (55) 
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FIG. 6. Transition rates of cooling A- (solid line) and heat- 
ing A+ (dashed line), for large cooperativity C 2> 1 and 
small cavity decay rate, k <C v. The arrows and circles mark 
the dressed states' frequencies and the dark resonances, in 
correspondence with Fig. [3] The parameters are 7 = 10^, 
k — v/lQ, g — 10V, 57l = 12v, Si — 5 opt = 47.5^, S C 2 = Si — v, 
ifi = 7r/3, (f>h = <j>c = 0. The transitions rates are given in 
units of a, Eq. (|57[) , and A in units of v. 

and it takes on the minimal value (m) s t,min ~ 1/C for 
Si = S op t and C_ <C 1. In this limit the cooling rate reads 
W ~ A_ = \e\ 2 Tf4:g 2 /-f. Figure [6] displays the transition 
rates A± for <5i = S opt . For this choice of the rate A- 
has a pronounced maximum at A = v. This maximum 
is displaced from the resonance around A w 2v by the 
the trap frequency ^, and corresponds to the case where 
the pump laser is resonant with the red sideband of a 
dressed state. The transition rate of heating is strongly 
suppressed at A = v due to three-photon resonance. 



C. Double dark resonances 

A peculiar characteristic of this specific setup is found 
when the mechanical effects are predominantly governed 
by the interaction with the cavity field. In this case the 
transition rates read 

A ± ^Aa ± ^ l2 \^ ± \ 2 + 2n\f^ ± \ 2 (56) 

with the amplitudes Tq^, Tq' , given in Eqs. (fl2"j) . 
These scattering amplitudes are connected to products 
of the functions T 1 and T K , Eqs. (|28|) . which determine 
the excitation spectra of the atom and cavity. 

Figure [7] displays the rates Ac,± for different values of 
the detuning Si and as a function of A. Cooling is found 
whenever Ac- > Ac+, that is, when the black curve 
overtops the gray curve. This is the case on the left 
edge of the broad resonance around A w Av. More in- 
teresting is the behavior around three-photon resonance 
A = Si — S C 2, marked by the black arrow. In particu- 
lar in the graphs b) and c), the rates Aq± oscillate with 
poles at a distance of the order of the trap frequency v, 
These oscillations are a consequence of Fano-like profiles 
in the scattering cross-section of the system and are vis- 
ible in the excitation spectra of atom and cavity. Here, 




A/u 

FIG. 7. Transition rates of cooling Ac,- (solid line) and heat- 
ing Ac,+ (dashed line) in units of a, Eq. (|57|) . as a function 
of A in units of v. The rates are evaluated from Eq. (|56(l for 
a) Si — <5 C 2/2, b) Si = 5 C 2 and c) Si — 3<5 C 2/2. The other pa- 
rameters are 7 = Wv, n = 2v, g — 20V, Q,l = 12z/, 5 C 2 = 20v, 
W2 = 1, if = it/3, <f>c = 0. The small arrows in the plots 
mark the positions of three-photon resonance (black) and of 
the corresponding sidebands, shifted by ±v (gray). 

these approximate double dark resonances lead to alter- 
nating cooling and heating regions around three-photon 
resonance. 

It is interesting to consider whether the parameter can 
be adjusted, to suppress both carrier and blue-sideband 
transitions, thereby increasing the cooling efficiency. An 
example for such an situation is the cooling scheme dis- 
cussed in Sec. IIVB 21 This could constitute a more ac- 
cessible realization of a cooling process, which so far has 
been proposed in a tripod level configuration (20| and for 
a two-level atomic transition confined in a resonator and 
at the node of the laser standing wave [34| . 

D. Large three-photon detuning 

The case where the control laser is far detuned from 
the atomic transition \gi) [e) corresponds to a closed 
two-level system consisting of the atomic levels | (72) and 
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|e) (in presence of a finite decay probability into state 
|gi) the control laser acts as rcpump). In our model, this 
case corresponds to taking large three-photon detunings 
5tp and large values of the detuning \Si\. In this limit the 
mechanical effects are due to coherent atom-cavity inter- 
action and to spontaneous decay, while the effect of the 
transverse laser coupling to the atom can be neglected. 
The atomic system thus reduces to an effective two-level 
transition coupled to a cavity which is pumped. This 
situation differs from the one considered in Ref. [HI, EH 
where the two-level atom is driven by a transverse laser. 
The case of a trapped two-level system whose center-of- 
mass motion is cooled in a driven cavity shows interesting 
features on its own, and its dynamics will be object of 
future investigations. 



V. DISCUSSION OF THE RESULTS 

In this section we report the graphs of cooling rate and 
the average vibrational number at steady state, (m) st . 
These quantities, defined in Eqs. (fM| and (f56")l . deter- 
mine the velocity of the cooling process and the final 
temperature, and are found by evaluating the rates in 
Eqs. ((57)) as a function of the physical parameters. For 
convenience, we report the cooling rate T scaled with the 
frequency 



a) 



-100 -50 50 100 



(57) 



The parameters we consider base in most cases on the 
experiment in Ref. (37| . Typical values for a, which are 



in accordance with Eq. (fl8|) . are of the order of 10 v to 

We first consider the three-photon resonance condition. 
This is achieved by appropriately setting the frequencies 
of the transverse laser and of the laser pumping the cav- 
ity. In this regime we expect that cooling is dominated by 
the mechanical effect stemming from the laser and cavity 
interaction, which mutually interfere, while the diffusion 
processes due to spontaneous emission are suppressed. 

Figure[8]displays the contour plots of the cooling rate T 
and of the corresponding stationary mean phonon num- 
ber (m) st as a function of the detuning A between probe 
and cavity, and of the detuning <5i between control laser 
and the atomic transition \gi) «-» |e). The solid line cor- 
responds to the condition given in Eq. (T4"Tj) . which relates 
A and 8\ and which maximizes the value of A- . Indeed, 
the lowest values of (m) st , and correspondingly the max- 
ima of the cooling rate, are found along the solid line. 
The largest value of the cooling rate T, Fig. [5] a), is 
found at the intersection between this curve and A = 0, 
namely, when the pump laser drives resonantly the cavity 
mode and the states \g2, 1) and \g±, 0} are hence also res- 
onantly coupled. This behavior can be understood con- 
sidering that the atom appears transparent to the cav- 
ity at three-photon resonance and in zero-order in the 
Lamb-Dickc expansion. In this regime, hence, the cavity 
resonance line, which scales the cooling rate, is given by 
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FIG. 8. a) Cooling rate V and b) mean phonon number (m) a t 
as a function of the detunings 5i and A when the three-photon 
resonance condition is satisfied, Stp = 0. The solid curves in 
a) and b) correspond to the resonance condition from Eq. (|47[) 
which maximizes A- (the dashed line reports the curve where 
A+ is maximum). In the hatched area the cooling rate is neg- 
ative and corresponds to the region where the atomic motion 
is heated, c) (m) st as a function of the detunings 5i for A 
fulfilling Eq. (g7j (solid line) and A = (broken line). The 
other parameters are 7 = 2u, n = I//4, g = 12^, J1l = 11.2v, 

if — 7r/3. 



a Lorentzian shape, and has its maximum when A = 0. 
Figure [8ji) displays the mean phonon number as a func- 
tion of 5± when A = (solid curve). This value is com- 
pared with the one found for A(d±), which satisfies Eq. 
(U7J and is indicated by the dashed curve. In both cases 
the system is driven at three-photon resonance. One ob- 
serves that the temperature is minimal for the values 
predicted by Eq. ([47|). 

Figure [5] displays cooling rate and mean vibrational 
number at steady state for the same parameters as in 
Fig. [8j but larger spontaneous emission rate 7 = 20f 
and smaller cavity decay rate k = O.OIza With respect 
to the results found in Fig. [SJ we observe that the cool- 
ing region about A = — v and b\ < becomes smaller. 
More striking is the appearance of a cooling region which 
stretches around A = v across all considered values of the 
detuning S± . This is understood in terms of the cooling 
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FIG. 9. Same as Fig. [8]but for k < i/. The plot c) shows cuts 
along the line A = (broken), A = v (dashed-dotted) and 
for A fulfilling Eq. (|47[) (solid line). The other parameters are 
7 = 20^, k = v/100, g = I2v, n L = 11.2v, if = tt/3. 



dynamics discussed in Sec. IIVB 21 In fact, for the con- 
sidered parameters the rate A_ is more than two orders 
of magnitudes larger than the rate A+. In this case, 
the cooling transition is enhanced because the pump is 
tuned to the red sideband of a resonance, while the heat- 
ing is suppressed because the corresponding transition 
coincides with a dark resonance. 

Figure [U) displays the mean phonon number as a func- 
tion of the detuning S\ for various values of A. The bro- 
ken line corresponds to the case A = and shows a simi- 
lar behavior as in Fig. [8l The thick dashed-dotted line is 
found for A = v: Here, not only diffusion is suppressed 
due to the three-photon resonance condition, but also the 
heating rate is strongly reduced. Indeed, the temperature 
is minimal for the value of Si fulfilling Eq. (|47|) at A = v. 
Further insight can be gained by inspecting the dressed 
states shown in Fig. [3] One finds that the cooling regions 
mostly correspond with the probe laser being set on res- 
onance with the red sideband of a narrow-line dressed 
state. We note, in particular, that the condition of opti- 
mal cooling can be identified in Fig. [3] at the intersection 
of the curve A = Si — 5 C 2 with the curve showing the fre- 
quency of the dressed state, shifted downwards by the 
trap frequency. This intersection determines the curves 
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FIG. 10. Mean phonon number (m) st as a function of the 
detunings A and 5i for fixed 8 C 2 = 20v. The solid curve 
corresponds to Re/(A + v) = 0, Eq. (|41[) . which maximizes 
A- (dashed curve reports Re/(A — v) — maximizing A+). 
The area in the box is magnified in the right plot. The other 
parameters are 7 = Wis, k, = 2v, g — 20v, J7l = 12v, S C 2 = 
20//, W2 = 1, <p = tt/3. 



of optimal cooling in Fig. [5] b) . 

We finally analyze the cooling dynamics when the 
three-photon resonance condition is not satisfied, re- 
stricting to some specific parameter regimes. Figure [TU] 
displays the mean phonon number at steady state as a 
function of A and Si at a fixed value S C 2 — 20V. The 
solid line corresponds here to the values of A and <5i ful- 
filling Eq. (|4ip . For comparison, the curve corresponding 
to three-photon resonance is reported and corresponds 
to the diagonal line. Minimal temperatures are mostly 
found along to this curve, demonstrating that EIT cool- 
ing is in most parameter regimes an efficient cooling 
scheme. Further parameter intervals where low temper- 
atures are encountered are close to the curves A = —30^ 
and A = 0. These regimes correspond to dynamics which 
can be understood in terms of Doppler cooling: the probe 
laser is here tuned on the red of the corresponding dressed 
state, whose linewidth is larger than the trap frequency. 
A strikingly different behavior is found close to the inter- 
section between the three-photon resonance curve and 
A = 0: this parameter region is reported in the inset. 
Here, we observe that heating and cooling regions al- 
ternate. The cooling dynamics is here determined by 
interference processes: The corresponding rates A± are 
displayed in Fig. 0(b), showing that the cooling rate is 
characterized by a double Fano-like profile. 



VI. CONCLUSIONS 

A cooling scheme has been presented and character- 
ized, which combines cavity quantum electrodynamics 
and electromagnetic induced transparency. Novel dy- 
namics have been identified, which may allow one to effi- 
ciently cool the atomic motion. These originate from the 
composite action of cavity, laser, and quantized atomic 
motion, giving rise to interference processes which allow 
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one to efficiently cool the motion. The cooling dynam- 
ics is robust against fluctuations of the parameters for 
most situations identified. Moreover, efficient ground- 
state cooling has been predicted for parameters which 
are consistent with the ones reported in Ref. (37j . 

The theory we presented here should also be applicable 
to current experimental setups, where the atom is con- 
fined in traps with different potentials for the relevant 
atomic states, as it is often the case for a dipolar trap 
setup. The fluctuations of the state-dependent trapping 
potential may be accounted for by a modified diffusion 
coefficient in the rate equations, as initial considerations 
suggest [llj]. Moreover, the theory can be extended to the 
case, when the motion is only strongly confined along the 
cavity axis. The perpendicular motion of the atom can be 
taken into account by averaging g over the corresponding 
slice of the cavity's mode function. 

Future work shall address the properties of the scat- 
tered light. For several cavity-based systems it has al- 
ready been demonstrated that photons can be emitted 
from a single atom in a controlled way [58l - l60j . Further- 
more, quantum interference sustains the creation of non- 
classical features in the emitted light [6l|, [62J. Therefore, 
the system here presented exhibits several requirements 
to be potentially operating as a quantum emitter. In 
general, given the large number of control parameters 
at hand, and the possibility to interface it with several 
fields, this system can result to be a very useful element 
for a quantum network [63} for which further building 
blocks based on elementary light-atom interaction have 
recently been realized [64| - [66j | . 
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Appendix A: Evaluation of the scattering rates 

1. Calculation of the transition amplitudes for 
cooling 

The transition amplitudes 

% = {<pt\ (v + V^-^V^j l^i). (Al) 

of a scattering process are evaluated using the resolvent 
of the Hamiltonian [o7| . The corresponding transition 



rates are 

11 f 

where the 5-function ensures energy conservation, and 
the sum goes over all relevant, different final states. The 
initial and final states and \ip{) have defined energies 
E{ and Ei with respect to the unperturbed Hamiltonian 
H , Eq. dlTJ). The total Hamiltonian U = H + V + O{if) 
contains the interaction part 

V = [Uq - H a ] +H!+W P 

= W 7 (0) + W K + Fx + F y } x + W P (A3) 

.7=1,2 

which includes the mechanical effects of laser and cavity 
F, Eq. (|2"5j) . and the mechanical effects denoted by F lj: 
Eq. (|26p of the interaction with the free radiation field in 
first order in r/. 

We regard scattering processes with initial state 

\<p i ) = \V st ,m;0 i J, (A4) 

where the atom and cavity mode are in the state | W st ) , 
Eq. (|23p . m phononic excitations arc present in the 
ccnter-of-mass motion of the atom, and the free radia- 
tion field is in the vacuum state. The initial energy is 
Ei = (m + \)hv + hA. We are interested in the pro- 
cesses which change the motional state of the atom by 
one phonon, corresponding to the final states, 

|^ ± ) = |* st ,m±l;lg) (A5a) 

|^' ± ) = |* st) m±l;lg). (A5b) 

These are states, where a photon is scattered into the 
modes of the electromagnetic field external to the res- 
onator, either by spontaneous decay along |e) f-> \g2) or 
cavity losses. 

If one neglects the decay along transition \gi) <h> |e), 
the states (|A5|) are the only final states for processes in 
first order of e and r\j where the motional state is changed 
by one quantum. We make the assumption 71 = to 
avoid multi-photon scattering at the transition driven by 
the strong control laser. This simplification allows for 
a less involved analytical treatment. Comparisons with 
numerical results show, that the results are not signifi- 
cantly altered in the parameter range we are interested 
in. 

The evaluation of the scattering amplitude 7fi, 
Eq. (|A1[) . using the initial and final states from Eqs. (|A4|) 
and (|A5p . is performed perturbatively in a Born series 
expansion 

% a (tpt\ [VGoiEJV + VGo(Ek)VG (Ei)V] \<pi) (A6) 
using the unperturbed resolvent 

Go(z) = 7^V' (A7) 
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containing the effective Hamiltonian 

Hf = H -ih~\e){e\-iHKa f a. (A8) 

which accounts for radiative losses. The evaluation of 
Eq. (|A6J) , using the initial and final states (|A4j) and (|A5|1 , 
reveals that only three processes described by the ampli- 
tudes 

l£ = (<p]' ± \F^G Wp\tp i ) (A9a) 
V* = (tf^lW^GoFGoWrltpi) (A9b) 
Tp'^ = (ip^IWk GoFGo W P \<pi) (A9c) 

contribute. Here, Go = Go (Si). Equation (|A9ap repre- 
sents the diffusive process due to the mechanical effect 
of spontaneous emission, whereas Eqs. ()A9b[) and (|A9cp 
stand for the lowest order mechanical effects of the con- 
trol laser and the cavity. 

With the specific definitions given in this appendix, it 
is easy to check that the scattering amplitudes, Eqs. (|A9[) . 
are the only contributions to light scattering including 
mechanical interaction up to first order in 77 and e. They 
are straightforwardly evaluated by substituting the vari- 
ous operators by the explicit expressions from Eqs. ([lip . 
(US), (ESI), ([26]) . The matrix elements of G (£i) are eval- 
uated in App. [Bj in particular Eqs. (|B1|) and (|B4j) with 
Ei = HA + h(m + \ One finds 

T± = h^T6^[e x ■ e^Jg^Mlkl), (A10) 

V !± = Kyfr+^g { £ ); %'H\k\), (All) 
= h^TdZg^f**^), (A12) 

after collecting all remaining factors - which are inde- 
pendent of the polarization e g and direction k/\k\ of 
the scattered photon and the vibrational quantum num- 
ber m — in the tilded quantities. The abbreviation S± 
yields 1 for the "+" case and otherwise vanishes. 

The corresponding rates are found using Eq. (|A2|) . 
where one has to sum over all relevant polarizations and 
wave vectors, 

R± = 2n(m + 5 ± )^ M [S x ■ %Js£ ) 'S(c\k\ ~ uj e2 )\f D \ 2 , 

k,e 

R^ ± ^2 7 :(m + 5 ± )^ {K \glll\ 2 S(c\k\~u Jc )\f^ ± \ 2 . 

k,e 

Moreover, we used wp rs luq ~ ^>ei 3> v to approximate 
the argument of the delta function, and w e 2 = uj e — LJ2 ■ 
The quantities 7d, 7p'' ± and Tp'^ in the latter equations 
are evaluated at the |fc| -value determined by the delta 



function. Taking into account the definitions [671 ] 

^ = \^9 { ff\ 2 5(hc\k\ - ^ e2 ), (A13) 

« = |El^fc?e| a W*|- fiw c), (A14) 

fc,e 

and A±(m + 5±) = i?f5 + #f * + ^f* 0] leads to ex- 
pression (|37| . where W2 was introduced as a geometri- 
cal factor characterizing the atomic dipole radiation pat- 
tern [23. 

2. Calculation of the excitation spectra 

The excitation spectra of the atom at rest give the 
rate of photon scattering as a function of the probe fre- 
quency wp either at the cavity output, giving S£ xc (uip), 
or from the resonance fluorescence of the atom integrated 
over the solid angle involving the modes external to the 
resonator, which we label by S^^cup). It has the form 
Si xc cx i? J (a;p), for j = k, 7, whereby i? J (a;p) is the 
rate of scattering a probe photon, assuming the atom- 
cavity system being in the initial state |<?2,0) at energy 
Ei = hup , and that the atomic motion is neglected. They 
are connected to the scatter amplitudes 

V = (tpi\W^Go{Ei)Wp\<pi). (A15) 

by Eq. (|A2|) . evaluated up to first order in the small pa- 
rameter e. These amplitudes correspond to processes, 
where, starting from the state Eq. (|A4j) . the cavity 
is excited by the probe laser as described by Wp. The 
subsequent evolution determined by Gq(E{) couples the 
states of the manifold {|e, 0), \gi, 0), \g2, 1)} due to the 
strong cavity-atom and control laser interaction. Gener- 
ally, all states of the manifold can be populated at this in- 
termediate state, and a photon can leave the atom-cavity 
system either by spontaneous emission (Wj = W J2 (0)) or 
by leaking through the cavity mirrors (Wj = W K ), lead- 
ing to the final state \ ff) = \g2,0,m; lg e )- We note that 
for spontaneous decay it is sufficient to take into account 
only W l2 (0), even for 71 7^ 0: Both scattering amplitudes 
only differ by a constant prefactor. 
Evaluating Eq. (|A15[) yields 

T< = ^g { £h{eAG { ol{hA)\g 2 , 1), (A16) 
TK = ^g^h(g 2 ,l\G^ t (hA)\g 2 , 1), (A17) 

where the required matrix elements of G^ p \(hA), cal- 
culated in App. [51 give the functions cosyJ r 7 (A), 
Eq. (|28bj) . and F K (A), Eq. (|28aj) . respectively. Using 
Eq. (|A2p and after summing over all relevant modes k, 
e of the emitted photon, one finds the excitation spectra 
in the form of Eqs. (|2"T1) . 



15 



Appendix B: Matrix elements of the resolvent 



subspaces of the state {|#2,0)} and the manifolds 



In this appendix the relevant matrix elements of the re- 
solvent Gq(z), Eq. (|A7[) . are calculated. We first rewrite 



Go{z) — G op t(z — H ext — -f/free) 

using G opt (z) = [z- i?opt] _1 and 



(Bl) 



^opt = tfopt - ih^\e)(e\ - ihnata. (B2) 



The operator G op t(z) has block diagonal form 

eff , 
Dpt J 



(P n H^ t P m = for n 7^ to) when considering the 



J 



M r 



{|fl 2j n),|e J n-l) > |5 1 ,n-l)} (n > 0) , (B3) 



where we denoted the corresponding projectors by Pq = 
92,0)(52,0| and P n , respectively. 

For the dynamics analyzed in this work it is sufficient 
to calculate the matrix elements belonging to the case 
n = 1. The relevant matrix elements are found by in- 
verting z — PiH£p t Pi, yielding 



mo 



|e,0)(e,0| {(<J c2 -Si + C)(« + C)} 



+152, l)( 9 2, 1| | (S c2 + i\ + C) (5 c2 - 6 X + C) - ^} 
+lffi,0)(5i,0| {(*«a + i| + c) (*« + - 5 2 cosV} 

+9 cos Lf — [|£H,0)(ff 2 ,l| +H.c] 
+(^ + C)^[|5i,0)(e,0|+H.c.] 
+gcostp(5 c2 -8i +C) [|g 2 ,l)(e,0| + H.c] 



(B4) 



All the matrix elements of Eq. (|B4[) are inversely pro- 
portional to the function 

/(c) = dct[c - PiH^p./h) = n cc - o 

f=±,o 

=(i« + C) { [<5c 2 + C + *|] («ca + C - *0 - X 



g 2 cos 2 ip{S c2 +C-S 1 ). 



(B5) 



Its poles are the complex eigenfrequencies Wj_ ff , a;® 5 of 



the non-Hermitian operator H°^ t , Eq. (|B2[) . reduced to 
the states within the manifold A^i. The real and imag- 
inary part of these eigenfrequencies are plotted in Fig. [3] 
and determine the relevant resonances of the scattering 
processes. 
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